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■ Abstract 

We study the infrared behaviour of the pure Yang-Mills correlators using 

■ relations that are well defined in the non-perturbative domain. These are the 
£2 . Slavnov- Taylor identity for three-gluon vertex and the Schwinger-Dyson equa- 
tion for ghost propagator in the Landau gauge. We also use several inputs from 

' lattice simulations. We show that lattice data are in serious conflict with a 

IT} • widely spread analytical relation between the gluon and ghost infrared critical 

exponeessencnts. We conjecture that this is explained by a singular behaviour of 
' the ghost-ghost-gluon vertex function in the infrared. We show that, anyhow, this 

discrepancy is not due to some lattice artefact since lattice Green functions sat- 
Qh. isfy the ghost propagator Schwinger-Dyson equation. We also report on a puzzle 

| concerning the infrared gluon propagator: lattice data seem to favor a constant 

non vanishing zero momentum gluon propagator, while the Slavnov- Taylor iden- 
• , tity (complemented with some regularity hypothesis of scalar functions) implies 

rN ' that it should diverge. 
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1 Introduction 



1.1 Generalities 

The whole set of correlation functions fully describes a Quantum Field Theory, as 
it is related to the S-matrix elements. In QCD or pure Yang-Mills theories Green 
functions are most often gauge dependent quantities which have no direct relationship 
with physical observables, the latter being necessarily gauge invariant. However, their 
indirect physical relevance is well known. In particular, long distance (or small mo- 
mentum) Green functions will hopefully shed some light on the deepest mysteries of 
QCD such as confinement, spontaneous chiral symmetry breaking, etc. In this paper 
we concentrate our efforts on the study of the gluon and ghost Green functions at 
small momentum in a pure Yang-Mills theory. Our tools will be Slavnov- Taylor (ST) 
identities, Schwinger-Dyson (SD) equations and also several inputs from lattice QCD. 

The infrared behaviour of Green functions has been extensively studied using differ- 
ent techniques, such as Schwinger-Dyson equations (see e.g. [3,5,24,25] and references 
therein), renormalization group methods [26], stochastic quantization [4,27]. These 
equations are exact consequences of QCD and can be easily derived using the path 
integral formalism. However their practical use reveals in most cases very difficult 
and one has to resort to a truncation which lessens the rigour of the method. One of 
the noticeable exceptions is the Schwinger-Dyson equation for the ghost propagator 
which contains only one integral and thus needs no truncation; this is the only one 
which we shall use in what follows. On the other hand, in order to exploit it in prac- 
tice, one usually has to make appropriate ansatze for the gluon propagator and the 
ghost-ghost-gluon vertex. 

We shall also use the Slavnov- Taylor identity which relates the three-gluons vertex 
to the ghost-ghost-gluon vertex in covariant gauges. Applying this relation in the non- 
perturbative domain will lead, under some assumptions about the infrared regularity 
of the dressing functions, to non trivial and surprising conclusions. 

Lattice simulations are, of course, another major tool to study small momentum 
Green functions. However this paper is not meant to be a standard "lattice paper" . We 
aim at using SD and ST to derive properties of the small momentum Green functions 
and we will use lattice simulations as valuable inputs in our theoretical discussion 
and as a check of some hypotheses. As we shall see, the outcome proves to be quite 
surprising and undermines some widely spread beliefs. 

In what follows we work in the Landau gauge, but some of the results we present 
are actually valid in any covariant gauge. Our notations are the following : 
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f abc r^ P (p,q,r) 
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rt c (p,k;q) = r bc (-ip u )g T^(p,k;q) (4) 

= f abc (-iPv)go 

■ [5^a(p, k; q) - q„k^b(p, k; q) + p v q,,c{p, k; q) + q u p^d(p, k; q) 
+ PvP^ e (P, k] q)] 

respectively for the ghost propagator, the gluon propagator, the three- gluons vertex 
and the ghost-ghost-gluon vertex 3 . All momenta are taken as entering. In eq. (@J) — p 
is the momentum of the outgoing ghost, k the momentum of the incoming one and 
q = —p — k the momentum of the gluon. F(p 2 ) and G(p 2 ) are the dressing functions 
of the ghost and gluon propagators respectively. We parameterise the propagators in 
the infrared by setting at leading order 



2\ 



9 \ a O 



G{p , , v 

fp 2 \ aF ^ 
F(p 2 ) — I — - ] , when p 2 is small, 

\tJ 

where A, 7] are some dimensional parameters. 

Let us make a brief and partial summary of the present predictions for otp^c ( 
[3-5,24-26,28,29]). We refer, for a very complete list of references, to the review by 
Alkofer and von Smekal [3]. All those references assume the relation 2otF + otG = and 
parameterise cxf and «g as 

a G = 2k sd . 

Different truncation schemes for the Schwinger-Dyson equations give ( [3,5,24-26, 
28,29]) 

k sd = 0.92 [29] 

k sd € [0.17,0.53] [5] , using a 2-loop perturbative input 
while another approach, ( [4,26]), predicts two possible solutions 

k sd = 1, or 
KsD = 0.59 

Lattice simulations give «g ~ 1. Note that «g = 1 corresponds to the gluon propagator 
being finite and non-zero at vanishing momentum; in other words, among the numbers 
we have just quoted, only the ones given in ref. [5] lead to a divergent infrared gluon 
propagator while the other values correspond to a vanishing one 4 . 



1.2 Numerical setup of the lattice simulations 

In this part we briefly describe the technical details of our numerical lattice simu- 
lations. 

3 We stick to the decomposition given in ref. [19] except for the arguments of the scalar functions, 
for which we keep the same order as in T itself 

4 After the completion of this paper our attention was drawn on ref. [6] which also predicts an 
IR-divergent gluon propagator and on refs. [7, 8] which lead to a finite non-zero one. 
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We use the standard Wilson action. For the SU(2) gauge group we have used 
lattices of size 32 4 and 48 4 with (3 = 2.3. This value of (3 corresponds to /?su(3) ~ 5.75. 
In the case of the SU(3) gauge group the simulation has been done on 32 4 , 24 4 and 
smaller lattices with (3 = 5.75 and (3 = 6.0. Those rather low values of j3 have been 
chosen because they allow measurements at small momenta. We have used periodical 
boundary conditions for the gauge field. The gluon propagator is defined as a mean 
value over gauge field configurations : 



The ghost propagator is calculated by the inversion of the discretised Faddeev-Popov 
operator (cf eq. l2*7j) . For this purpose we have used the conjugate gradient algorithm 
with the source 



where V is the number of lattice points. The role of the 1/V terms is to eliminate the 
zero modes of the Faddeev-Popov operator (corresponding to global gauge transforma- 
tions) in order to allow its inversion in the orthogonal subspace (cf ref. [2]). All lattice 
data have been extrapolated to the continuum as described in ref. ( [2]). A detailed 
report on all numerical results is presented in the same reference. 

2 Constraints on ap and ac from the ghost 
Schwinger - Dyson equation 

There is a widely used relation between otp and Qg, referred to in the following as 
R Q , which comes from the scaling analysis of the Schwinger-Dyson equation for the 
ghost propagator ( [4]) and states that in four-dimensional space one has 2aF + ctG = 0. 

We have attempted to test this relation on lattices with the characteristics indicated 
above. We plot in fig.[T]the quantity F 2 (p)G(p) as a function of p. If the relation R a was 
true, this quantity should be constant in the infrared domain. One can see that it is 
not the case: F 2 G goes to zero at small momenta. On the other hand the ultraviolet 
(UV) behaviour is exactly the expected one. The same trend is already visible at /3's 
larger than ours : in refs. ( [12,13]) it is mentioned that F 2 G might decrease as the 
momentum approaches zero 5 at (3 = 6.0 and (3 = 6.4. The same authors have even 
reported on the same effect in the unquenched case [11]. Very recently Ilgenfritz et 
al. have published in ref. [16] results which go in the same direction (although the 
conclusions they draw thereof differ from ours). 

It was suggested in [13] that Gribov copies might induce significant changes in the 
infrared behaviour of the ghost propagator. Could this explain our findings for F 2 G ? 
In ref. [17] the accurate lattice gauge fixing (choosing the "best copy", corresponding 
to the lowest value of ||A^ 9 ' ) ||) seems to lessen the infrared divergence of the ghost 
propagator [14, 15, 17], implying a further increase of the drop of F 2 G in the infrared, 
while the gluon propagator is known to be only slightly affected by the presence of 

5 In ref. [12] it is assumed that the vertex function stays constant in the zero momentum limit, 
in which case F 2 G is proportional to the strong coupling constant a s in the MOM scheme based on 
ghost-ghost-gluon vertex. 



< A«(x)A b M > . 
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lattice Gribov copies. Furthermore, we shall see that the SD equation is closely related 
to a lattice-SD equation which is a mathematical identity, valid independently of the 
choice of the Gribov copy. Thus the possible influence of Gribov copies on propagators 
cannot explain the behaviour of F 2 G at small momenta. 
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Figure 1: F 2 G from lattice simulation for the SU{3) (left, 32 4 , /? = 5.75) and SU(2) 
(right, 32 4 and 4:8 4 ,Psu(2) = 2.3) gauge groups. The /3's are chosen so as to give the 
same lattice spacings : 1.2 GeV" 1 . If the relation 2a f + «g = was true this quantity 
should be constant in the infrared domain. One clearly does not see this behaviour. 

The rest of this section is aimed at understanding this disagreement between lattice 
simulations and the theoretical claim (R Q ). We will first revisit the proof of the latter 
in order to identify all the hypotheses needed and submit each of them to a critical 
analysis. In the next-to-following section we will discuss a special writing of the ghost 
Schwinger-Dyson equation, in a form which involves only Green functions instead of 
vertices and can thus be directly tested on the lattice. 

2.1 Revisiting the relation between ap and olq 

We will examine to what extent the proof of R Q is compelling, using the Schwinger- 
Dyson equation for the bare ghost propagator which can be written diagrammatically 

as 



V 



o 

k 



V 




i.e. 
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(^ (2) )^) = (?) 

-9lf acd fe b f I 7^^eH?)(^)r^(-g,^g-A:)^)(G( 2 ))2(g-A;), 

where we use the notations of eqs. (013) and go is the bare coupling constant. 
This integral equation is written in terms of bare Green functions. It can be cast 
into a renormalized form by multiplying (resp. F^ 2 \T) by Z^ 1 (resp. Z% ,Z\) 

and $q by Z^ 2 = and multiplying the k 2 term by Z%. The integral therein is 

ultraviolet divergent but one can check that the cut-off dependence is matched by the 
cut-off dependence of Z g and of the Z3 factor multiplying k 2 . Later on we will only use 
subtracted SD equations such that the UV divergence is cancelled as well as the Z^k 2 
term. These subtracted SD equations hold both in terms of bare and renormalized 
Green functions without any explicit renormalization factor. 

Let us now consider (J7|) at small momenta k. The ghost-gluon vertex may be 
expressed as 

qv'^u'u{-q, k;q-k) = quli^q, k) + (q - k) u H 2 {q, k) (8) 
where, using the decomposition (JJJ), one gets: 



Hi(q,k) = a(-q,k;q - k) - (q 2 - q-k) (b(-q,k;q - k) 
+ d{— q, k;q — k)) + q 2 e(-q, k;q — k)) ~ 

~ a(— q, k;q — k) — q 2 (b(—q, k;q — k) + d(—q, k;q — k) — e(—q, k;q — k)) 
H 2 {q,k) = {q 2 - q-k)b(-q,k;q - k) - q 2 c(-q,k;q - k) (9) 

In the Landau gauge, because of the transversality condition, H 2 does not con- 
tribute. Thus, dividing both sides of eq. (JJJ) by k 2 and omitting colour indices, one 
obtains 



1 _, , 2 N f d*q F(q 2 )G((q-kY 
F(k) L + 9 ^ C J (2vr)4 I q^q-kf 



q 

(q-k) 2 



Hfak)). (10) 



2.2 What does the "non-renormalization theorem" exactly 
say? 

A widely used statement, known as the "non-renormalization theorem", claims that, 
in the Landau gauge, the renormalization constant Z\ of the ghost gluon vertex is 
exactly one. Note that there is no reference to a particular renormalisation scheme. 
Formulated in this way, this claim is wrong. Let us first state and then explain below 
what is true in our opinion : 

1) There is a true and very clear statement which can be extracted from Taylor's 
paper (the argument is given below), ref. [1]. 

rp^ e (-p,0;p) = -z/^ (11) 
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i.e. there is no radiative correction in this particular momentum configuration (with 
zero momentum of the ingoing ghost) 

2) This entails that Y abc ' Bare {p, k; q) is finite whatever the external momenta, and 

— -MS ~ MOMj 

that therefore Z\ = 1. In addition, we get also trivially Z\ =1, where MOM/, 
refers to the configuration of momenta in equation In general, in other schemes, 
there is a finite renormalisation, and this is why we do not adopt the misleading ex- 
pression " non-renormalization theorem" . 

3) In particular, one finds in the very extensive calculations of radiative corrections 
at least two cases of MOM schemes where there is a finite renormalisation (and cer- 
tainly many more) : MOM s in the notations of ref. [19], and the symmetric MOM 
scheme. For the latter, we give the proof below. 

The essence of Taylor's argument is actually very simple. In a kinematical situation 
where the incoming ghost momentum is zero, consider any perturbative contribution 
to the ghost-gluon vertex. Following the ghost line in the direction of the flow, the 
first vertex will be proportional to the outgoing ghost momentum p M , i.e. to the 
gluon momentum — p^. In the Landau gauge this contribution will thus give upon 
contraction with the gluon propagator D^p). Therefore the only contribution to 
remain is the tree-level one. In other words the bare ghost-gluon vertex is shown to be 

equal to its tree-level value in these kinematics : Y abc,Bare (— p, 0; p) = —if abc Pfj,. This 
result has been checked by means of a direct evaluation to three loops in perturbation 
theory by Chetyrkin. In our notations : 



H^p, 0) + H 2 (p, 0) = 1. 

Note that in the Schwinger-Dyson equation EH only Hi is present, and the theorem 
of Taylor does not tell that Hi(p, 0) = 1, as seems assumed in many Schwinger-Dyson 
calculations, where it plays a crucial role in the proof of H a . 







o 
o 



6- 







Figure 2: The kinematical situations considered below. The left diagram (0- momentum 
incoming ghost) corresponds to Th below which is known to be equal to one. The right 
one (0-momentum gluon) corresponds to T g and leads to a non-trivial p 2 - dependence 

As an illustration of our point 3), let us quote the formulas from the appendix 
of ref. [19], reduced to the situation we are interested in = 0,nf = 0). The two 
dressing functions T h (resp. T g ) are defined by r^f c (— p, 0,p) = —if abc T h (p) (resp. 

r™ fec (— p,p, 0) = —if abc T g (p) ) and correspond to the kinematical situations depicted in 
the left (resp. right) part of fig. (J2J). We have already mentionned that Th is exactly 
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one, but this does not hold for T g and, indeed, one has at three loops : 
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It is then easy to find the p 2 -dependence : 



r„ 



-pMS 

g 
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In ref. [28] the non -renormalization theorem is understood as the statement that 
the vertex reduces to its tree-level form at all symmetric-momenta points in a symmet- 
ric subtraction scheme. However this statement is not supported by a direct evaluation. 
Using the one-loop results of Davydychev (ref. [20]) one gets in a symmetric configu- 
ration the value 



r?>, k; g)| p2=fcW= , 2 = -if abc (l + ^ ^(9 + ^)) + q, g ^(3 + ^)}(13) 

with = ^c/ 2 (f ), cud = 1M9 ■■■■ 

According to ref. [28] the coefficient of p M should be one. The presence of a s in 
the above formulas implies on the contrary that the vertex will in general depend on 
the momenta : using the results given in the appendices of ref. [19] one finds for the 

leading ^-dependence -if ahc {^{%f log(g) ((9 + §0)^ + (3 + |0)g M )}. This de- 
pendence is logarithmic, as is expected in a perturbative approach. Furthermore, in 
ref. ( [28]) it is supposed that the vertex function takes the form (q 2 Y(k 2 ) m ((q — k) 2 ) n 
with the restriction I + m + n = 0. One should note that this last condition corre- 
sponds in our notations to ctr = (cf. section ()2.3j) below). This restriction comes 
from the assumption that the symmetric vertex is equal to 1 for any p 2 , which, as we 
have just seen, is actually not the case. Therefore we shall adopt a more general point 
of view and keep open the possibility of a non perturbative effect leading to a singular 
or vanishing limit of H\ when q — > 0. We should mention that the problem of the 
p 2 -dependence of the ghost-gluon vertex has already been addressed in refs. [33,34]. 
However these authors work under the condition R a which appears not to be satisfied 
by our lattice data. 



2.3 A subtracted Schwinger-Dyson equation 

Let us now consider two infrared scales ki = k and k 2 = Kk. Calculating the difference 
of eq. (fTUj) taken at scales k x and k 2 and supposing for the moment that dp ^ one 
obtains 



1 



F(k) F(nk) 



OC (1 — K 



-2a ? 



){k 



2\-a f 



9oN c 



d 4 q / F(q 2 



(27T) 



k 2 q 



X 



G((q - kf)H x {q, k) G((q - Kk^H^q, nk) 



(14) 



{{q-kf 



{{q-Kkf 
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We now make the hypothesis that there exists a scale go such that 



G(q 2 ) ~ (q 2 ) 



2\a G 



F(q 2 



(g 2 )^, for g 2 <g 2 . 



Similarly, we suppose that Hi can be written for q 2 , k 2 < q 2 as 

q ■ k k 2 * 



Hi(q,k) ~ (q 2 r 



q 2 ' q 2 



or as 



H\(q, k) 



((q-k) 2 r h 2 



q ■ k k 2 
g z g z 



where the scalar functions h 12 are supposed to be regular enough (i.e. free of singu- 
larities worse than logarithmic) and expandable in Taylor series for k — > They are 
obviously invariant under any simultaneous rescaling of both q and k. The exponent 
ar gives the leading critical behaviour of H\ on q. 

Thus we rewrite (I14j) by rescaling — > A A; with A chosen so that (AA;) 2 <C g 2 an d 
splitting the integral in the r.h.s. into two parts 



f A 




l#<4 (2^) 4 





[ d% 









(15) 



In Ji, since (AA;) 2 -C g 2 , we can substitute the infrared approximations (0) for G and 
F . Ii is infrared convergent if : 

o>f + ar > — 2 IR convergence at g 2 = 

«G + «r > —1 IR convergence at (g — A;) 2 = and (g — kA;) 2 = (16) 

We shall suppose in the following that these conditions are verified. We then obtain, 
performing the change of variable q — > Xq and writing generically h for h 12 : 



-Tl(A) ~ / \ 2 ( Q F+ a G+or) 



d A q 



2\ap+ar- 



-i({k-qf 



k' 2 



4 w 

((? - k)T G - 2 h (f , J) - ((g - Kk) 2 ) aG ~ 2 h (k% k 2 %) 



(17) 



The point we have to keep in mind is the fact that the upper bound of the integral 
goes to infinity when A — ■> 0. This potentially induces a dependence on A whose 
interplay with the behaviour explicitly shown in (fT7|) we must check. In this limit, the 
convergence of the integral depends on the asymptotic behaviour of the whole integrand 
for large q. In particular, the leading contribution of the square bracket in eq. (fT7|) 
behaves as 



-,2\a a -2 



q 



2« G -6 



because the terms in q ■ k, being odd under g M — > — g M , give a null contribution under 
the angular integration in eq. (fT7|). Thus, assuming the conditions (fTB^) are satisfied 
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CclSG 


ap ^= u 
ap + «g + ar < 1 


ap ^= u 
ap + «g + ar > 1 


ap — 
ac + ar < 1 


ap — 
«g + ar > 1 


IR 

v^onvei g,ence 

Conditions 


cxp -\- ar ^> — 2 
«G + «r > —1 


ck ^ — (— ar > — 2 

«G + Otr > —1 


ar .> — z 
«G + ar > — 1 


ar -> — £ 
ac + ar > —1 


constraints 


2a f + «g + «r = 


ap = —1 


excluded 


none 



Table 1: Summary of the various cases regarding the a coefficients 



the integral Ii(A) is guaranteed to be convergent when q — > (or k — > g) and its 
asymptotics for small A is given by 



h (A) 



/•go /a 

^2(a G +a F +a r ) / 2(a P +a! G +a: r )-3 ^ ^2(a G +a F +a r ) 

Jo 

if «g + op + ar < 1 

A 2 / dq g2(a F+aG +ar)-3 _ A 2 if « G + «p + Q<r > 1 
JO 



because in both cases the integral on the momentum q is finite and does not depend 
on A in the limit A —>■ 0. 

Let us now consider I2. Its dependence on A is explicit in the factor 

G {{q - Xk) 2 ) H{q, \k) _G{{q- Xnk) 2 ) H{q, Xnk) 

({q-\kyf ((g-A^) 2 ) 2 



which stems from the substitution k — > A k in (JHJ). Clearly, this quantity can only be 
even in A : any odd power of A would imply an odd power oiq-k whose angular integral 
is zero. Since the integrand is identically zero at A = and the integral is ultraviolet 
convergent, it is proportional to A 2 (unless some accidental cancellation forces it to 
behave as a higher even power of A). 

So, if the first of the conditions (JTHJ) is verified, it follows h + 1 2 « A 2(aF+ ° G+Qr) , 
else I\ + I2 ~ A 2 . Comparing this to the left hand side of eq. (fTU|) 

(19) 



F(Xk) F(\Kk) 
we then conclude : 

«i? + «G + «r <1 ==> 2ai? + ac + ar = 

ap + ac + ar >1 =>■ ap = —1 (20) 

In the particular case ap = the leading term in the l.h.s. of eq. ()14)) is identically 
zero. We are left with the subleading one which, pursuing the same argumentation, we 
suppose to be proportional to k 2 . Then the argument is the same as in the previous 
case except for the power of A in the r.h.s. of eq. (fTTJj) which becomes equal to 2. It 
results that the case «g + ot < 1 is now excluded while the case «g + or > 1 provides 
no extra constraint. 

The various possibilities which have appeared in this discussion are summarised in 
table From this table it appears that only the triple condition that 



ap 7^ 0, ar = 0, ap + ac < 1 
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does actually imply the standard statement that 2ap + «g = . However, the plot 
in FigQ indicates a behaviour 2a p + > 0, indicating that at least one of these 
conditions is not fulfilled. 

Let us assume for the moment that aj? + etc > 1 and ctr = 0. Then 2aF + ac > 0, 
in agreement with fig^ an d etc > 2. However the possibility that «g be greater than 2 
is unambiguously excluded by the lattice simulations so that the hypothesis has to be 
rejected. Furthermore, as we shall see, one can derive from the Slavnov- Taylor identity 
relating the three-gluons and ghost-ghost-gluon vertices the inequality etc < 1 if one 
assumes that some of the scalar form factors of these vertices are regular when one 
momentum goes to zero. 

We now consider the hypothesis 2a f + etc + «r = with «r < to comply 
with the lattice indications of figffl This implies that some of the scalar factors of 
the ghost-ghost-gluon vertex are singular in the infrared. We shall turn back to this 
possibility in the concluding remarks of this section (|2.4.4J) . The question is whether a 
non-perturbative effect could generate a non vanishing ar- 

A direct lattice estimate of the ghost-ghost-gluon vertex would be welcome. How- 
ever this is a difficult task. A direct measurement implying a zero momentum ghost 
is impossible since the corresponding Green functions are singular because of the zero 
modes of the Faddeev-Popov operator. A careful limiting procedure implying very 
small external momenta has to be performed. This study is under way. In between we 
propose a simpler check based on another writing of the SD equation (JZj) in terms of 
a pinched Green function which can be directly checked on the lattice. This will also 
allow to control the lattice artefacts and to rule out the hypothesis that the problem 
encountered in fig. [T] is simply due to a lattice artefact. 

2.4 Green function formulation of the ghost SD equation 

We will now rewrite the ghost SD equation using only propagators and the ghost- 
ghost-gluon Green function. In this form, its validity can be tested on the lattice, 
because what one directly calculates in lattice simulations are Green functions, not 
vertex functions. If we consider the loop integral in Eq. (JZJ) we can see that it is 
nothing else but a ghost-ghost-gluon Green function in which the left ghost leg has 
been cut and where the gluon and right-hand ghost have been pinched onto the same 
point in configuration space. We shall see that this quantity is directly accessible from 
lattice data without making any specific assumption about the behaviour of the vertex 
function. 

The interest of this approach is that it will help us to throw a closer look at the 
compatibility of the lattice simulations and the SD equations. Indeed, as we have 
just seen, we are facing a contradiction between, on the one hand the lattice estimate 
of F 2 G (fig. and, on the other hand, the relation R a which is derived from the 
SD equation complemented by a regularity assumption (ar = 0) suggested by 
perturbation theory. Therefore we feel the need to directly confront lattice calculations 
with SD. The form of SD which is presented in this subsection allows such a direct 
confrontation and, this form being closely related to a lattice SD equation which is just 
a mathematical identity, we are in a good position to trace back any discrepancy. 

In the next subsection we present the continuum limit derivation of the Green 
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function formulation of the ghost SD, as well as its lattice derivation. 

2.4.1 Continuum limit case 

For a given gauge field configuration A^ = A®t a , the Faddeev- Popov operator in 
covariant gauges is given by 

M(x,y) lom i=(d lt D li )5(x-y) = (A + ig Q d^) 5\x - y) = F^~\a, x, y), (21) 

and it is equal to the inverse of the ghost correlator in the background of the gluon 
field A, ^iconf- The subscript means here that the equation is valid for any given gauge 
configuration. This can be written 

- y ) = M lconf (x, z)F^ mt (A, z, y), (22) 

where a summation on z is understood. Expanding M according to (|21|) we get 

8{x -y) = A(x, z)FgUz, y) + ig 8® (a»(x)f£Ua, x, y)) , (23) 

valid for any gauge field configuration. Performing the path integral one gets the mean 
value on gauge configurations 

S(x -y)= A(x, z) (F^Uz, y)> + i^W (A^F^A, x, y)). (24) 

Of course {F^ nf (z, y)) is nothing else but the ghost propagator defined in equation pf. 

The averages (F^(x,y)} and djj C \A tl (x)F^(x, y)) are invariant under translations 
so that one can replace the derivative d^f 1 —>■ —0$ . We take x = and perform the 
Fourier transform on the y variable (note that there is no tilde on -4^(0) ) 

1 = -p 2 F\p) - g p,(A,(0)F^ o JA,p)). (25) 

Finally we get 

F(P 2 ) = 1 + 9oj^f ab %Al(0)F^(A,p)). (26) 

One caveat is in order here. Eq. (J2J) implies an ultraviolet divergent integral which is 
matched by renormalization constants. This divergence is of course also present in eq. 
(124)) via the local product of operators at x. In section l2~T1 this divergence was canceled 
by subtracting two terms as is apparent in eq. ()14|) . In the following this divergence 
is regularised by the lattice cut-off. To perform the connection with the discussion in 
section 12.11 it will be necessary to perform an analogous subtraction when using the 
form (J2HJ). 

2.4.2 Lattice case 

We now repeat the same steps as in the preceding paragraph for the lattice version 
of the Faddeev-Popov operator 

K = F^jU^y)' 1 = [ Sa » - <W) - S*(x - e „) {5 yn - 5 y , x ) + 

+ \r hc [^(z)<w, - a;(x - e^^J ' , 

(27) 
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where 

S;\x) = -±Tr [{t a ,t»} (U,{x) + Ul(x))} 

A (rr\ - U^(.x)-Ul(x) _ 1 T U^(x)-Ul(x) W 
' n H\- L J 2 N 2 ' 

in which U^(x) denotes a standard link variable 6 , and e M is a unitary vector in direction 
/i. We define 

= (S« b (x) (5 x , y - 5 y , x+e J - Sf(x - e,) (5 yn - S y , x )) . (29) 

The appearance of as the appropriate discretisation of the usual Laplacian operator 
A is dictated by the gauge invariance of the original Yang-Mills action, which imposes 
that the standard V operator be replaced by its covariant version and by the specific 
form - 9foTr(^ links U 9 ) - of the functional to be minimized in order to fix the Landau 
gauge. 

Then, multiplying M^ b y by F^(x,y) from the right, one obtains : 

N 2_i TrA u(v> Z ) F SL( U ;*,«) = s v,u- 
c fabc , . (30) 



2(iV c 2 - 1) 



^) F iconf (U; y + e„ u) - A%y - e,)F^{U- y - e„ it) 



This is an exact mathematical identity for each gauge configuration which must actually 

(2) 

be fulfilled by our lattice data since our -F lcc , nf are computed by means of an explicit 
inversion of the Faddeev- Popov operator. From this fact results an important feature 
which we wish to stress : since eq. valid in any configuration, its consequences are 
free of any ambiguity originating from the presence of Gribov copies. Upon averaging 
over the configurations one gets 



^ F —[T r (A u ( y ,z)F^ n Az,u)) = 6 y , u 

jabc 



(2) 

(31) 



C 

£abc 

( 2 ) ba {TT „, , „ „.\ At (a, „ \ T?( 2 ) ba 



2( y2 } L (i;(|/)C(^ V + e„ u) - A^y - e^F£Z(U, y - e„ «)) 

Of course, this averaging procedure depends on the way chosen to treat Gribov's prob- 
lem : the particular set of configurations over which it is performed depends on the 
prescription which is adopted (choice of any local minimum of A 2 , restriction to the 
fundamental modular region...) and, consequently, the Green functions will vary but 
they will in any case satisfy the above equation. Like in the continuum case, after 
setting y to zero, a Fourier transformation with respect to u gives: 



^Tr j p - U (&u(0, z)F£ljU, z, «)) 



£abc _ (32) 

Note that although eqs. (jHOJ) and (jSH) have to be exactly verified by lattice data 
eq. (|3*2*|) does only approximately (within statistical errors) since it relies on transla- 
tional invariance, which could be guaranteed only if we used an infinite number of 
configurations. 

6 Note that the definition of given in 1|28J) differs from the naive one by a factor igo. This is the 
reason for the presence of i and the absence of go in eq. I|31|> as compared to eq. l|lUfl • 
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Eq. fK->2l) is a discretised version of (J2EJ)- Therefore any lattice correlator, satisfying 
(J32J) . should also satisfy ()26|) up to non-zero lattice spacing effects. Among the various 
sources for such effects the use of the specific Au discretisation of the Laplacian opera- 
tor in the l.h.s. deserves some comments. The gauge fields present in the S^ b (x) terms 
in eq. (ffilj) generate in Ajj the so-called "tadpole" diagrams such as in fig. El According 



2 

a 



Figure 3: Example of the terms in the Schwinger- Dyson equation on the lattice. 

to the philosophy developed by Lepage and Mackenzie in ref. ( [23]) the tadpole contri- 
bution can be estimated by a mean field method. Using the average plaquette P (for 
(5 = 6.0 P ~ 0.5937) one predicts a tadpole correction factor oc P - ^/ 4 ) ~ 1.14. These 
terms disappear in the continuum limit but they do so only very slowly : the tadpole 
corrections (1 - plaquette) vanish only as an inverse logarithm with the lattice spacing. 
This is to be contrasted with the corrections arising in the r.h.s which are expected to 
be of order a 2 . 

2.4.3 Checking the validity of the SD equation on the lattice 

Since, as we have just mentionned, eqs. (j3ll|) and (}3Tj) are mathematical identities, 
there is in principle little -if any- to learn from a verification of eq. (J32|) . except for 
the verification that our configurations are actually in the Landau gauge. On the other 
hand one thing we wish to be reassured about is the possible role of lattice artefacts 
in the discrepancies we have noticed. 

We begin with a comparison of the continuum r.h.s of eq. ()26|) with the l.h.s. of (}3*2"|) . 
Both sides are plotted in fig. (JH). 

The agreement is impressive. Should we on the contrary use eq. ()26j) itself we 
observe a clear disagreement between the two sides of the equation. What we thus 
learn is that the major part of the discretization artefact comes from Ajj. In fact, 
our lattice data show that Au ~ A/1.16 almost independently of the momentum, see 
fig. El This is in good agreement with the correction factor of 1.14 obtained from 
Lepage-Mackenzie's mechanism [23]. To conclude the tadpole effect explains almost all 
of the discrepancy observed when trying to verify ()26)1 . One can also understand why 
this discretisation artefact is so large : this is due to the slow logarithmic vanishing of 
the tadpole corrections. 

It results from this discussion that the lattice artefacts cannot be blamed for the 
violation of the relation R a observed in fig. : the tadpole effect has been seen to 
produce a corrective factor almost constant in p, thus unable to explain an error in the 
power behaviour. 

2.4.4 Concluding remarks 
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Figure 4: Check of the validity of the SD equation on the lattice. 
Left : ^{Au&Tr PgUp)) (circles) is plotted vs. 1 +g O ^ i f abe (A%(0)F^(A, p)> 

(squares). Right : ^(A^Tr F( 2 )(p)> -ffo^/^^^ffl^p)) is compared 
to 1. 
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Figure 5: Lattice computation of A/Ajj in Fourier space as a function of the momen- 
tum. The statistics is poor for technical reasons. Note also that the region above 
7r/2 (~ 3 GeV) is affected by strong discretisation effects. 

We want to emphasize at this stage that our lattice data both satisfy the properly 
discretised SD equation and violate the relation R a . The most likely way out we can 
think of is that the hypothesis = is not verified, i.e. that Hi is singular when all 
momenta are small . One more argument in favour of this explanation is provided by 
a direct numerical examination which shows that, in the eventuality of neglecting the 
momentum dependence of Hi, the SD equation is satisfied in the UV but badly violated 
in the IR. The results of this comparison are given in fig. El- The details of the method 
can be found in the appendix. 

The possibility of a singular behaviour of Hi in the infrared has already been 
considered by various authors, for example in ref. [21] in the framework of exact renor- 
malization group equations. It is known from perturbation theory ( [1]), that the sum 
Hi + if 2 (cf eq. 0) is exactly equal to 1 at k — (this has been checked explicitly 
in perturbation theory up to fourth order in ref. [19]). This is the argument which is 
usually called for when advocating «r = 0. However it does not prevent the contribu- 
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Figure 6: Testing the truncated ghost SD equation on the lattice. The data correspond 



to a 577(3) simulation at (3 — 6.4 and V = 32 



3.66 GeV). 



tions of the scalar function b to H% and H 2 , which cancel out in the sum, from being 
singular in this limit. Actually, from dimensional considerations, one concludes that 
b must be of dimension —2. At k = the only dimensional quantity involved (at the 
perturbative level) is q, which means b(q, 0; —q) oc 1/q 2 . This singularity is removed 
by the kinematical factor in front of b in Hi and if 2, but this would no longer be the 
case for k 7^ if one had more generally b(q, k; —q — k) oc 1/q 2 . In any case our results 
appear to plead in favour of a divergent ghost- ghost- gluon vertex in the infrared domain. 



Slavnov- Taylor identity and the infrared behaviour 
of the gluon propagator 



Another non-perturbative relation that can be exploited is the Slavnov- Taylor iden- 
tity. It can be used to constrain (under some hypotheses) the infrared exponent for the 
gluon dressing function. In the preceding section we have explored the consequences 
of the very strong assumption that Hi is regular when all its arguments go to zero and 
we have shown that this assumption is not tenable when the lattice data are taken into 
account. We shall now make the weaker hypothesis that the scalar factors present in 
the decomposition (J3J) of the ghost-ghost-gluon vertex are regular when one of their 
arguments go to zero while the others are kept finite and exploit the Slavnov- Taylor 
identity under this assumption to derive constraints ona F and a<> 
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The Slavnov- Taylor ( [1]) identity for the three-gluon function reads 



p X T XfML ,(p,q,r) = 53 P J (hvr 2 - r x r v )T \p(r,p; q)- 



G(r 2 ] 
F(p 2 ) 
G{q 2 ) 



(33) 



(<W<? 2 - q\q^\u(q,p; r). 



One then takes the limit r — > while keeping q and p finite. The tensor structure of 
rV/(p, k; q) has been recalled in eq. (jlj). We adopt the following notations for particular 
kinematic configurations: 



a>z{p 2 ) = a(-p,p;0) 

ai (p 2 ) = o (0, &i(p 2 ) = 6(0, di(p 2 ) = d(0, -p;p). 



(34) 



In the present case of one zero momentum the three-gluon vertex may be parameterised 
in the general way as ( [19]) 



Fftvpip, -p, 0) = (25 ^pp - S^pPv - SpvPfj) Ti(p 2 )- ( 8^ - 



PyPv 
p2 



P P T 2 (p 2 ) +Pp,PvP P T 3 (p 2 ) . 

(35) 



Now, exhibiting the dominant part of each term of eq. ()33|) we obtain: 

Ti(q 2 )(q^ - q 2 5p U ) + q 2 T^(q 2 )q^q v + T]i^(q,r) = 
(ai(g 2 ) + r 1 (q,r))(8p, l/ r 2 - ryv)l + + {b 1 (q 2 ) + r 2 (q,r))q tl (r 2 q u - q.rr u ) 



G(r 2 ) 



+ {bi{q 2 ) + di(q 2 ) + r 3 (q, r))r M (r 2 ^ - q.rr u ) 



+ 



I F {{g+rf) 
+ G(9») 

where r lj2 ,3 and r/i i2 verify 



(36) 



limr^g, r) = lining, r) = limr 3 (g, r) = 

i — >0 r— »0 r^O 

lim?7i^(g,r) = lim 772^ (g, r) = 



(37) 



Identifying the leading terms of the scalar factors multiplying the tensors q^q u and 
QfiQv — q 2 8p U we obtain the usual relations ( [19]): 



f) = ifW? 2 ) 



Ti{t 

n(q 2 ) = 0. 



(38) 



Using these relations we see that Eq. ()36j) implies: 



. F{p 2 



lim 



G(r 2 



a x (q )(r 5^ - Tpr v ) + bi(q )(r g M g„ - r.qq^) 







(39) 



Thus one sees that if ai(g 2 ) 7^ or b\ 7^ (and, indeed, one knows from perturbation 
theory that at large momenta a\ = 1, cf. [1, 19]) (J33|) can only be compatible with the 
parameterisation (j5J) if 

a G < 1. (40) 
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We can also, instead of letting r — > 0, take the limit p — > of Eq. ((331) as i s done 
in [19]. The dominant part of the l.h.s. of (|33J) is : 

{25^p.q - p^q u - Puq^Tiiq 2 ) - [8^ - { ^)p.qT 2 {q 2 ) + p.qq^vT z {q 2 ) 

The r.h.s. is the product of F(p 2 ) with an expression of at least first order in p. T\ 
and T 2 being different from zero we can conclude in this case that ap < 0. 

Let us repeat here that all these considerations are valid only if all scalar fac- 
tors of the ghost-ghost-gluon and three-gluons vertices are regular functions when one 
momentum goes to zero while the others remain finite. Under those hypotheses one 
obtains important constraints on the gluon and ghost propagators - namely that they 
are divergent in the zero momentum limit. 

3.1 Lattice results 

The results for a G and ap from lattice simulations are presented in figs. (JTJ) and (JJJ). 
We present in TableElthe values of the coefficients for a fit of the form 7 (q 2 ) a {\ + ^q 2 )- 



Group 


Volume 


P 


a G 


ap 


SU(2) 


48 4 


2.3 


1.004 ±0.015 


-0.087 ±0.015 


SU(2) 


32 4 


2.3 


0.968 ±0.011 


-0.109 ±0.014 


SU(3) 


32 4 


5.75 


0.864 ±0.016 


-0.153 ±0.022 



Table 2: Summary of the fitting parameters for the F and G functions 

The fits have been performed without using the point p 2 = even in the case of the 
gluon propagator where it is known. Its inclusion would have forced a G to be equal to 
1. In any case one sees on fig. [7|that the point at p 2 = available in the S'[/(2),32 4 
case is compatible with the fit. 

For SU(2) and the larger volume the value obtained for a G is compatible with 1. 
The situation is less clear for SU(3), but in this case data with the larger volume (48 4 ) 
are lacking. Moreover we have to take into account our experience from previous studies 
of the gluon propagator where we have always observed that the gluon propagator goes 
continuously to a finite limit in the infrared region. A very detailed study of the gluon 
dressing function and specially of its volume dependence at k = has been performed 
by Bonnet et al. (cf ref. [30]). This study shows that a value a G = 1 is compatible with 
the data (the dressing function shows no signal of discontinuity in the neighbourhood 
of zero) and that no pathology shows up as the volume goes to infinity. 

We conclude that lattice data seem to contradict Zwanziger's result (G( 2 )(0) = 0), 
[22] and most probably also the predictions derived from the Slavnov- Taylor identity 
((j( 2 )(0) infinite). Like in the preceding section a possible way out of this contradiction 
could consist in dropping the regularity assumptions which have been made in the 
course of the proof. 

7 A term of the form /iq 2 is clearly needed in order to describe a situation like the one in fig. J3 
left) where G^(p 2 ) seems to go to a finite limit when p goes to zero 
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Figure 7: G (2) (p 2 ) from lattice simulation for SU(2) (left) and SU(3) (right) . /3 SU(2) = 
2.3 and /3 SU (3) = 5.75. the volumes are 32 4 and48 4 for SU{2) and 32 4 for SU(3). 

4 Discussion and Conclusions 



4.1 Discussion of the validity of SD and ST equations in the 
infrared 

Schwinger-Dyson equations and Slavnov- Taylor identities are valid non-perturba- 
tively. However some care is needed mainly due to Gribov ambiguities in the gauge 
fixing procedure. One needs a well defined QCD partition function for a given gauge 
(in the following we concentrate on Landau gauge). Many different non-perturbative 
prescriptions for the quantisation of non-abelian theories have been suggested : in- 
tegration only over the absolute minima of the J A 2 functional (Gribov fundamental 
modular region) [22], summing of copies with signed Faddeev- Popov determinant [9], 
stochastic quantisation [31], etc. All these prescriptions correspond to different valid 
gauge fixing procedures. 

In lattice numerical simulations two main methods are of practical use: The algo- 
rithm which minimises the functional is stopped at its first solution, which is a local 
minimum, or one takes the smallest local minimum among a given number of trials on 
the same gauge orbit. One is sure to be inside the Gribov region, never to be inside 
the fundamental modular region. 

The question is whether SD and ST relations are valid in these gauge fixing schemes. 
It is argued in ( [4]) that the Schwinger - Dyson equations are valid under different 
quantisation prescriptions provided that the Faddeev-Popov determinant vanishes on 
the boundary of the integration domain. However, the partition functions will differ 
and hence, the Green functions will be in general different solutions of the the SD 
equations. 

The Slavnov- Taylor identities may be derived from the QCD partition function 
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using the gauge invariance of the action and stating that the gauge fixed path integral 
is invariant after a change of variables which corresponds to a gauge transformation. 
It is then clear that the Slavnov- Taylor identities remain valid whichever gauge fixing 
procedure has been followed. However this proof has to be taken with care in the 
presence of singularities of the partition function or of Green functions. 

One more general comment is in order. As we have already stated, no Green 
function with a vanishing ghost momentum can be defined on the lattice since the zero 
modes of the Faddeev-Popov operator are discarded. For the same reason no source 
term for the zero mode ghost field is allowed. More generally we do not know of any 
non perturbative way to define such Green functions. We cannot prove that this means 
a divergence of the Green functions with one ghost momentum going to zero, but we 
can suspect that it is the case as lattice indicates for the ghost propagator 8 . Indeed 
the close-to-zero modes of the Faddeev-Popov operator are strongly influenced by the 
Gribov horizon and, for sure, very different from any perturbative result. This casts 
also some doubt about the use of ST identities or SD equations in the case of a zero 
ghost momentum. 

4.2 Conclusions 

We have tried in this work to put together various inputs in order to clarify our 
understanding of the infrared behaviour of the pure Yang-Mills Green functions. Our 
findings can be summarised as follows : 

1. The lattice Green functions contradict the common lore according to which 2a F + 
a G = 0. The present situation is that 2a F +a G > 0, i.e., the product F 2 (k 2 ) G(k 2 ) 
tends to for k — > 0. From what we observe concerning the evolution of the curves 
with the size of the lattice, it is difficult to imagine how a further increase of the 
volume could eventually revert this tendency. 

2. The result 2a F + a G = (R a ) which is contradicted by lattice data, is usually 
claimed to be derived from the ghost SD equation. Our results seem to cast 
some doubt on the validity of these derivations, based on the assumption of a 
trivial ghost-ghost-gluon vertex ("naive approximation"). Indeed, we do verify 
that the properly discretised SD ghost equation on the lattice is well satisfied. We 
conclude that the lattice data seem to prove that this "naive approximation" is 
invalid, and that there exists in the ghost-ghost-gluon vertex a non-perturbative 
infrared singularity of the form (k 2 ) ar which has been neglected in the standard 
analysis and which leads to the replacement of R a by 2a F + a G + C(r = (ccr < 0). 

3. Regarding the ghost correlator it results from its very definition that it cannot be 
defined at zero momentum. Its lattice values at low momentum appear to be in 
favour of a divergent dressing function (a F < 0), as is also suggested by the above 
theoretical arguments. However the divergence is much slower (a F e [—-15, — .1]) 
than what is usually reported. This is in agreement with the conclusions one can 
draw from the Slavnov Taylor identity. 

4. In relation with points 2) and 3) above it is worth insisting on the fact that the 
Schwinger -Dyson equation by itself is not sufficient to determine the behaviour of 

8 Notice that we cannot claim anything about the behaviour of vertex functions when one ghost 
momentum goes to zero. 
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the ghost-propagator and of the ghost-ghost-gluon vertex. Different treatments 
of the Gribov copies lead to different infrared solutions (se also ref [17] in this 
respect), all of which fulfill the SD equation. 

5. As for the gluon propagator the situation is much less clear. Three sources of 
information are available and give contradictory results. 

- The Slavnov Taylor identity, supplemented by regularity assumptions for 
the ghost-ghost-gluon vertex functions, points towards a divergent infrared 
behaviour. 

- Our lattice data indicate a finite limit when the momentum goes to zero. 
This trend is very clear in the SU(2) case although less compelling for SU(3). 
The fit has been performed by excluding the point at p = but the latter, 
when known, is compatible with the extrapolated value. These results agree 
with the ones of our previous studies as well as with the findings of the other 
lattice groups who have studied this matter ( [30]) which include the point 
p = and impose therefore = 1. 

- Zwanziger's result [22] states that the gluon propagator vanishes at k=0 but 
a fully satisfactory proof of a continuous vanishing as k — > is still lacking. 

We are unable for the moment to settle this point in a totally unambiguous way. 

This set of conclusions raises some questions. First of all further studies are still 
under way in order to fix the issues related to point 4). As to point 2), a direct lattice 
study of the ghost-ghost-gluon vertex at low momenta is desirable, although difficult. 
Such a study has recently been performed for SU (2) in a specific kinematical situation 
(zero- momentum gluon) ( [32]). The precise relevance of this special case to the points 
we have considered remains to be clarified. 

Note added 

After the completion of this work we realised that the particular situation (ap = 0) 
which we have mentionned in the fourth column of table ^ but not fully investigated 
might provide a good agreement with the lattice data while complying with the con- 
straints stemming from the SD equation. This possibility is discussed in a further 
publication ( [35]) . Let us stress that this solution too is compatible with the non- 
vanishing of 2a p + etc 
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6 Appendix: Testing the naive approximation of 
the ghost SD equation 

The simplest approximation of the ghost SD equation (fTT)|) corresponds to the case 
H 1 (q,k) = lVq,k: 

g 2 N c f d\ (F{q 2 )G{{k-qf){k-qf-kW 



1 + 



(41) 



F(k) k 2 J (2tt) 4 V q 2 (k-q) 2 (q - kf 

Strictly speaking this equation, written in this way, is meaningless since it involves UV 
divergent quantities but a corresponding meaningful renormalized version can be given 
(see the caveats about eq (JZJ) in subsection (|2.1jl ). We want to check whether lattice 
propagators satisfy it. According to perturbation theory, it should be true at large k. 
Lattice propagators are discrete functions, and thus one has to handle the problem of 
the numerical evaluation of the loop integral I in ([41)1 . Let us express the integrand in 
terms of q 2 and (k — q) 2 



d*q F 2 (q)G(k — q)i(k — q) 2 + (k 2 ) 2 + (q 2 ) 2 - 2k 2 q 2 q 2 + k 2 



(42) 



(2tt) 4 q 2 (k-q) 2 I 4 4(Jfe - q) 

Then we write 

I = h + 1 2 + h + h + h + h, 

each Jj corresponds to one term in (|42|). 
All these integrals have the form 

Ii = Ci(k) J ^-JM)K{k-q). 

The convolution in the r.h.s. is just the Fourier transform of the product at the same 
point in configuration space: 

d 4 q 



fi(q)hi{k -q) = F+ F_(/ 4 )NF_(/k)N (*) 



where F_ (/)(#) is an inverse and F + (f)(k) a direct Fourier transform. Thus, in order 
to calculate the integral I from discrete lattice propagators one proceeds as follows: 

1. calculate {fi}(p) and {hi}(p) as functions of F(p), G(p),p 2 for all i 

2. apply the inverse Fourier transform F_ to all these functions and get fi(x) and 
hi(x) 

3. compute the product at the same point fi(x) ■ hi(x) 

4. apply the direct Fourier transform F + to fi(x) ■ hi(x) 

The calculation of Fourier transforms involves a Hankel transformation which is nu- 
merically evaluated by means of a Riemann sum 

f(r) = (2n)- 2 \\r\\- 2 £ Mrp t )p 2 M±lhA {pt _ pi _^ po = . (43) 
i=i 

The inverse transformation is done in the similar way. In practice, because of the 
lattice artefacts which become important at large p the summation has to be restricted 
to p < p max ~ 2.2 instead of the "ideal" value 2n. 
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Errors There are three important sources of errors : statistical Monte-Carlo errors 
for F(q 2 ) and G(q 2 ), the bias due to the integral discretization and the truncation 
of the p-summation to values lower than p max . The second one is dominated by the 
neglected contribution coming from the UV cut-off of the integral (the integration is 
performed on some ball -8(0, L) instead of IR 4 ). Let us estimate the error on the Fourier 
transform of the product in such a case: 

F + (f(x)g(x))(k) = [ d 4 pd*q6 e (k-p-q)f(p)g(q), (44) 
Jb(o,l) 

where the e— approximation to delta function is 5 e (p) = f B , Q ^ d A xe^ x,p \ Considering 
e small enough we can integrate on q around the point (k — p), obtaining finally: 

(/ * 9)l « (/ * 9)oo + Vol(fl(0, e(L))) • (/ * V^oo 

This gives us an estimation of the error coming from the UV cut-off. As for the last 
source of error, it may be neglected because of the following argument : the integral 
is logarithmically divergent, therefore the neglected behaves as log (|^) — log (^f^) = 
lofff-^-). Thus it remains finite as a goes to zero and gets smaller and smaller as 

v pmax 7 

compared to the part actually computed. 

Results: We still have to face the same problem we have already encountered in 
section (|2. 4.3(1 . namely that the lattice Faddeev-Popov operator involves the non trivial 
discretisation Ajj ofj4ie Laplacian operator. This is taken into account by means of 
the substitution of Ajj(p 2 )/p 2 to the "1" term in the l.h.s of equation (|4*T|) We present 
on (Fig. [BJ) the result of the numerical integration described above. We have chosen 
for this purpose the data set from the simulation with the gauge group SU (3) at 
(3 = 6.4, V = 32 4 . One sees that the equality is achieved at large momenta, but in the 
infrared the naive approximation of the ghost SD equation fails. 
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